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The Baetis vernus group (Ephemeroptera: Baetidae) 
of northernmost Europe: an evidently diverse 
but poorly understood group of mayflies 


Eino Savolainen, Marcus K. Drotz, Per-Ola Hoffsten & Anssi Saura 


1. Introduction 


Savolainen, E., Drotz, M. K., Hoffsten, P.-O. & Saura, A. 2007: The Baetis ver- 
nus group (Ephemeroptera: Baetidae) of northernmost Europe: an evidently di- 
verse but poorly understood group of mayflies. — Entomol. Fennica 18: 160- 
167. 


The Baetis vernus group is known to be a taxonomically difficult group of may- 
flies. Morphological characters such as the mandibles or gills of nymphs overlap 
widely. We report here that enzyme electromorph frequencies also fail to show 
reproductive isolation in many cases. Baetis macani has been divided in the liter- 
ature into two forms, a lentic one with broader gills living in stagnant waters, 
while the lotic one with narrower gills lives in flowing water. We show that the 
two evidently represent reproductively isolated taxa, B. macani Kimmins = B. 
bundyae auct. nec, while the one with broader gills appears to be a new species. 
Baetis macani, B. subalpinus, B. liebenauae and B. vernus form a cluster that is 
genetically not divisible according to the recognized species boundaries. Nymph 
gill length differentiates, however, B. macani from B. subalpinus and B. vernus. 
We call for studies on population structuring and molecular taxonomy to resolve 
the taxonomy and ecology of this important group of mayflies. 
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biomonitoring of environmental quality (e.g. 
Fialkowski et al. 2003). 


The genus Baetis is geographically widespread 
and many species are abundant, in particular in 
northernmost Europe. Because the adult stage is 
short, nymphs are the most easily collected and 
biologically the most interesting stage in the life 
cycle. Baetis nymphs are now extensively used in 


Linnaeus (Linné 1767, p. 901) recognized the 
difficulties inherent in assigning the short-lived 
mayfly adults to species. The adults and nymphs 
of the northern European representatives of the 
vernus group (Miller-Liebenau 1969) of Baetis 
are difficult indeed to separate morphologically. 
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For a broader overview on the relationships 
within Baetis the reader is referred to Jacob 
(2003). Engblom (1996) gives an authoritative 
key to the species but when one is examining 
field-caught material there is much geographical 
and ecological variation in characters. To give an 
example, Bauernfeind and Humpesch (2001) re- 
viewed the central European fauna. They note 
that the gills of B. macani, vernus and tracheatus 
nymphs are considerably longer when they live in 
lake and pond habitats in comparison with flow- 
ing water. 

The vernus group Miiller-Liebenau (1969) 
comprises the following north European species: 

Baetis vernus Curtis 1834 exhibits much mor- 
phological variation (Engblom 1996). It has a 
wide transpalearctic distribution. The nymphs 
live in flowing water (Miiller-Liebenau 1969, Ja- 
cob 2003). 

Baetis subalpinus Bengtsson 1917 is a north- 
em species that has not been found outside of 
Fennoscandia. In northern Sweden the nymphs 
live in water flowing with intermediate velocity 
(Miiller-Liebenau 1969, Jacob 2003). 

B. macani Kimmins 1957 was originally de- 
scribed from Finnish Lapland from an arctic/al- 
pine pond (Kimmins 1957). Macan (1957) de- 
scribed the nymphs from the same pond and, in 
addition, from two large lakes from central Fin- 
land and central Sweden. The nymphs live both in 
lentic and lotic habitats, typically at outflows 
from lakes (Jacob 2003). 

Baetis bundyae auct. nec is a northern Euro- 
pean form with narrow gills. Engblom (1996) 
gives small ponds above the tree line as the typi- 
cal habitat in Europe. Baetis bundyae Lehmkuhl 
1973 was described from tundra ponds in arctic 
Canada. It has also been recorded in flowing wa- 
ter elsewhere in arctic Canada (e.g. Cobb & 
Flannagan 1980, Harper & Harper 1997) and has 
a distribution extending south along the Rocky 
Mountains. Morihara and McCafferty (1979) 
designated the Canadian Baetis bundyae 
Lehmkuhl as a subspecies of the European B. 
macani. The subspecific status is based on geo- 
graphical isolation, North America versus Eu- 
rope; the morphological variation is subject to 
“overlapping character states between the two 
subspecies”. There is certainly much overlap in- 
side Europe and ecological confusion as well. To 


give an example, Degerman ef al. (1994) give 
separate distribution maps for B. bundyae, B. ma- 
cani and B. macani-bundyae in Sweden. 

Baetis liebenauae Keffermiiller 1974 was de- 
scribed from the river Odra in western Poland 
(Keffermiiller 1974). In Fennoscandia it was first 
recorded from southernmost Sweden in 1979 
(Engblom 1996) and from northern Finland by 
Savolainen and Pulkkinen (1987). It has a distri- 
bution covering almost all of Europe and Asia 
Minor. All records come from flowing water (Ja- 
cob 2003). 

Baetis tracheatus Keffermiiller & Machel 
1967 was described from Poland (Keffermiiller & 
Machel 1967). Savolainen and Saura (1996) re- 
ported it from northern Finland. Both in Poland 
and elsewhere in Europe the nymphs live in lakes 
and slowly flowing water among water plants 
(Miller-Liebenau 1969, Jacob 2003). 

In addition, we have used nymphs of Baetis 
rhodani Pictet 1843-1845, B. niger (Linnaeus 
1761) and Acentrella lapponica (Baetis lapponi- 
cus) Bengtsson 1912 as an outgroup. Within the 
genus Baetis, B. rhodani belongs to the subgenus 
Rhodobaetis and B. niger to the subgenus Nigro- 
baetis (Jacob 2003). 

The existing keys for larvae are not always 
useful. We use here enzyme electrophoresis and 
allele frequencies as a criterion of reproductive 
isolation. Using an extensive material covering 
several years, we aim at seeing whether the 
nymph forms found in northernmost Europe re- 
present bona species. Specifically, we wish to 
find out whether the lotic and lentic forms as- 
signed to a single species are conspecific and re- 
present phenotypic responses to different envi- 
ronmental conditions or whether there are one or 
several ecologically isolated species and whether 
the nymph forms assigned to different species in a 
single habitat represent reproductively isolated 
entities. 


2. Material and methods 


The author Savolainen collected mayfly nymphs 
from the localities shown in Fig. 1 in northern 
Finland. Table 1 gives the names of the species, 
the names of the localities, the days when col- 
lected and the numbers of individuals assayed by 
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electrophoresis. In addition, the author Hoffsten 
collected nymphs from a locality from which 
Engblom had collected on 5.VII.1980 and had 
identified as Baetis macani with broad gills at 
Lake Bastuträsk, Västerbotten, Sweden (this is 
sample number 50 in the Table 1). The nymphs 
were kept in deep freeze until assayed for starch 
gel electrophoresis. The nymphs were assigned to 
species following the key of Engblom (1996), 
which relies on mandibles in particular. Nymph 
mandibles are subjected to wear, a circumstance 
that makes species identification difficult (see the 
discussion section for details). 

The electrophoretic methods are described by 
Savolainen et al. (1993). The following enzymes 
gave repeatable results: acid phosphatase (Acph), 
esterases (Est), glutamate oxaloacetate trans- 
aminase (Got), Peptidase (Pep), phosphoglucose 
isomerase (Pgi) and Triosephosphate isomerase 
(Tpi). In addition, all systems mentioned by Sa- 
volainen et al. (1993) were tried but the remain- 
der was regarded as insufficiently repeatable or 
they were monomorphic in all species studied 
here. Samples from different populations were 
run on the same gel so that the samples served as 
controls to each other. 

To get an overview of the extensive temporal 
and spatial variation disclosed with enzyme elec- 
trophoresis, one needs a method to visualize the 
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Fig. 1. Sampling localities of mayflies 
in northernmost Europe. 


results obtained from different gene pools. The 
neighbour-joining algorithm of Saitou and Nei 
(1987) is particularly suited for analyzing this 
kind of large data set. It does not assume constant 
rates of evolution and approximates the minimum 
evolution method (Page & Holmes 1998). The 
unbiased genetic distance of Nei (1978) was cho- 
sen as clustering algorithm. Phenogram construc- 
tion was done using Version 1.0 of the GDA pro- 
gram (Lewis & Zaykin 1996). We have only in- 
cluded populations that gave repeatable results 
over the set of enzyme systems in constructing 
the phenogram. 

We revise the species names in the discussion 
section. To avoid confusion, we use the revised 
names in the results. 


3. Results 


Distinguishing Baetis species on morphological 
criteria has been notoriously difficult. This diffi- 
culty is also reflected in the genetic differentia- 
tion: one can not recognize clear-cut species 
boundaries among B. subalpinus, B. vernus, B. 
liebenauae, B. tracheatus and B. macani (Fig. 2). 
The lentic form of Baetis macani with narrow 
gills had only moderate but consistent frequency 
differences against the lotic form with broad gills. 
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Table 1. Sample number, Baetis species name, site number (water course, parish, coordinates according to the 


uniform grid system), habitat (B = brook, R = river, L = lake), collecting date, and number of nymphs sampled. 
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macani 
macani 
macani 
macani 
macani 
macani 
macani 
macani 
macani 
macani 
macani 
lapponica 
lapponica 
liebenauae 
liebenauae 
sp. nov. 
sp. nov. 
sp. nov. 
sp. nov. 
sp. nov. 
sp. nov. 
sp. nov. 
sp. nov. 
sp. nov. 
niger 
rhodani 
rhodani 
rhodani 
subalpinus 
subalpinus 
subalpinus 
subalpinus 
subalpinus 
subalpinus 
subalpinus 
subalpinus 
subalpinus 
subalpinus 
subalpinus 
tracheatus 
tracheatus 
tracheatus 
tracheatus 
vernus 
vernus 
vernus 
vernus 
vernus 
vernus 
vernus 
undescribed 


1 Paskalomaoja, Savukoski, 7532:549 
1 Paskalomaoja, Savukoski, 7532:549 
2 Sarakko-oja, Savukoski, 7529:551 


3 Parakainen, headwater, Savukoski, 7535:550 


4 Parakainen, Savukoski, 7534:549 
5 Ahdaskuru, Enontekiö, 7678:251 


3 Parakainen, headwater, Savukoski, 7535:550 


6 Pälkevaara, Kittilä, 7539:389 

8 Pauttiselkä, Sodankylä, 7558:510 
9 Lutto, Inari, 7593:518 

10 Kaunispäänoja, Inari, 7595:516 
5 Ahdaskuru, Enontekiö, 7678:251 
12 Myllyjoki, Kittilä, 7527:407 


11 Niemisjoki, | rapids, Kuhmo, 7115:650 


12 Myllyjoki, Kittilä, 7527:407 

13 Luosujärvi, Kolari, 7495:376 

14 Riikonjärvi, Kittilä, 7519:417 

15 Ylläsjärvi, Kolari, 7495:385 

16 Särkijärvi, Muonio, 7538:371 

17 Törmäslommol, Muonio, 7543:374 
18 Rukajärvi, Kuusamo, 7340:599 
19 Äkäslompolo, Kolari, 7505:378 
20 Kallavesi, Kuopio, 6977:531 

22 Kelujärvi, Sodankylä, 7487:501 
11 Niemisjoki, Kuhmo, 7115:650 
23 Karigasniemi, Utsjoki, 7702:463 
11 Niemisjoki, Kuhmo, 7115:650 

7 Hillajänkänoja, Kolari, 7503:373 
12 Myllyjoki, Kittilä, 7527:407 

24 Kutuniva, Muonio, 7543:374 

25 Markkajoki, Enontekiö, 7598:379 
11 Niemisjoki, Kuhmo, 7115:650 
11 Niemisjoki, Kuhmo, 7115:650 
26 Keimiö-oja, Kittilä, 7540:380 

12 Myllyjoki, Kittilä, 7527:407 

27 Torpanoja, Kolari, 7495:385 

21 Jeesiöjoki, Kittilä, 7508:438 

28 Ylisenoja, Sodankylä, 7572:509 
29 Kiilo-oja, Sodankylä, 7583:517 
21 Jeesiöjoki, Kittilä, 7508:438 

21 Jeesiöjoki, Kittilä, 7508:438 

21 Jeesiöjoki, Kittilä, 7508:438 

21 Jeesiöjoki, Kittilä, 7508:438 

30 Orankijoki, Pello, 7427:373 

21 Jeesiöjoki, Kittilä, 7508:438 

21 Jeesiöjoki, Kittilä, 7508:438 

21 Jeesiöjoki, Kittilä, 7508:438 

21 Jeesiöjoki, Kittilä, 7508:438 

31 Mannajoki, Kolari, 7489:357 

21 Jeesiöjoki, Kittilä, 7508:438 


32 Bastuträsk, Skellefteä, 719900:169900 


Tr DUIDNIVMDNIVNVNUVVDNVDVVVVDVVUNUITTTITITTITT TT VDVOVDUUDTUUUTUUUUU 


5 August 1996 
5 August 1997 
5 August 1997 
5 August 1997 
5 August 1997 
8 August 1997 
8 August 1998 
10 August 1998 
22 July 1999 
23 July 1999 
24 July 1999 
8 August 1997 
5 August 1995 
7 August 1998 
10 August 1998 
3 August 1995 
4 August 1995 
5 August 1995 
6 August 1995 
6 August 1995 
2 August 1996 
7 August 1997 
6 August 1998 
27 July 1999 
7 August 1998 
3 July 1997 

7 August 1998 
21 July 1999 
5 August 1995 
6 August 1995 
7 August 1995 
3 July 1997 

7 August 1998 
10 August 1998 
10 August 1998 
10 August 1998 
22 July 1999 
27 July 1999 
27 July 1999 
2 August 1995 
3 August 1996 
6 August 1997 
9 August 1998 
2 August 1995 
2 August 1995 
3 August 1996 
6 August 1997 
9 August 1998 
21 July 1999 
22 July 1999 
25 July 1998 
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In contrast to the overall quite diffuse allelic dif- 
ferentiation seen in other Baetis taxa, the form 
with broad gills has a clearly circumscribed gene 
pool. Therefore we designate it here as sp. nov. 
Accordingly, the two macani forms represent re- 
productively isolated taxa. The allele frequency 
data are available in the form of an electronic ap- 
pendix (Drotz et al. 2006). 

Nevertheless, within the B. vernus group one 
could observe certain clear-cut genetic differ- 
ences. The Myllyjoki population number 29 mor- 
phologically identified to be B. subalpinus was 
fixed for a rare allele at Tpi that differed from the 


33 subalpinus 
44 vernus 
35 subalpinus 
4 macani 
14 liebenauae 
15 liebenauae 
40 tracheatus 
3 macani 
36 subalpinus 
39 subalpinus 
8 macani 
38 subalpinus 
7 macani 
48 vernus 
29 subalpinus 
1 macani 
2 macani 
6 macani 
11 macani 
9 macani 


Fig. 2. A neighbour-join- 
ing phenogram based 
on the clustering algo- 
rithm of Nei (1978) 
showing genetic dis- 
tances among the pop- 
165.105 ulations studied, based 
: on the allelic variation at 
the enzyme loci studied 
(see the electronic ap- 
pendix, Drotz et al. 
2006 for allele freguen- 
cies). Note that some 
populations are ex- 
cluded from the 


10 macani 
24 sp. nov.* 
32 subalpinus 
31 subalpinus 
42 tracheatus 
47 vernus 


18 sp. nov. 
21 sp. nov.* 
23 sp. nov.* 
22 sp. nov.* 
43 tracheatus 
20 sp. nov.* 
19 sp. nov.* 
17 sp. nov.* 


45 vernus phenogram, mainly be- 
5 macani cause of small sample 
28 rhodani size or the lack of data 
26 rhodani for a certain enzyme sy- 


stem. * undescribed 
species, see text. 


12 A. /apponica 
13 A /apponica 


one seen in other samples from the same locality 
(12). The locus is in general monomorphic within 
a species and different alleles fixed at nearby pop- 
ulations indicate reproductive isolation. The 
Bastuträsk sample from central Sweden (site 32) 
also differed from the northern Finnish ones far 
more than the Kallavesi population from central 
Finland did, even though geographically it was 
closer to the northern samples. Baetis tracheatus 
had allele freguencies indistinguishable from sp. 
nov. in one case. 

Outside of the B. vernus group one locus or 
two loci were often almost or completely fixed 
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for species-specific alleles that, accordingly, can 
be used as species diagnostic characters. Acen- 
trella lapponica differs from the others at the 
Acph and Pgi loci. It has the allele 85 fixed in Pgi, 
while all other species but B. rhodani have very 
low freguencies of this allele. The sample of 
B. niger was small indeed, but the three individu- 
als had species-specific alleles fixed at Est-2, Pep 
and Tpi loci. Baetis rhodani differs from the 
others at Acph and Pgi. 


4. Discussion 


4.1. A complex without distinct 
species boundaries? 


Baetis vernus group is known to be a morphologi- 
cally difficult cluster of species (Mäller-Liebenau 
1969, Jacob 2003). This difficulty is also seen in 
the allele freguency data generated by enzyme 
electrophoresis (Drotz et al. 2006). We may take 
as an example the large temporal variation in B. 
tracheatus from the same locality (21) sampled in 
1995, 1997 and 1998. One sample (43) is very 
close to the sp. nov., while two others (40 and 42) 
cluster with B. macani and B. liebenauae. In con- 
trast, the samples of sp. nov. come from an exten- 
sive area. 

Baetis macani, B. subalpinus and B. vernus 
cluster with each other. They all come from flow- 
ing water. The electrophoretic methods that we 
have used do not show consistent allele frequency 
differences that would follow the slight morpho- 
logical differences that characterize these taxa. 
Morphologically B. macani has consistently long 
gills and is a rather coherent taxon in comparison 
to B. subalpinus and B. vernus, which show char- 
acter overlap (cf. Bauernfeind & Humpesch 
2001, Engblom, pers. comm.). Species discrimi- 
nation is based mainly on mandibles, even though 
Miiller-Liebenau (1969) had pointed out that 
characters such as surface of tergites, epicranial 
sutures, bristles of femora etc. can be useful. We 
have observed that patterns of mandibular wear 
vary among individual nymphs, this again may be 
site specific. Miiller-Liebenau (pers. comm.) has 
stressed that because of wear, one should not look 
at the mandibles in use but at the next, unused pair 
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that is visible under the nymphal cuticle. 

The small genetic differentiation between B. 
macani/subalpinus/vernus that we have observed 
may be a consequence of introgression, i.e gene 
flow across species boundaries. We have studied 
only nymphs. The short-lived adults are, how- 
ever, the agents of reproductive isolation. Baetis 
adults swarm in July—August over low vegetation 
in a vertical fashion from 0.5 up to 4 m from the 
ground from afternoon until evening. The males 
will try to mate with about everything that flies in 
the vicinity (Savolainen 1978). They swarm close 
to the water, where the nymphs live, so that 
swarming either next to flowing or stagnant water 
results in ethological isolation. Savolainen et al. 
(1993) have shown that swarming behaviour is 
associated with genetic differentiation in may- 
flies. 

On the other hand, local selection such as wa- 
ter velocity (Bauernfeind & Humpesch 2001) 
could be expected to modify the flexible nymph 
morphology. Molecular studies using both 
mtDNA and nuclear markers would help solving 
these issues. Williams et al. (2006) have, using 
the CO1 gene of mitochondrial DNA found that 
Baetis rhodani is made up of a complex of cryptic 
species. In our opinion their results strengthen the 
conclusions that we have reached here: there is 
much ecological and geographic differentiation 
going on among Baetis maytlies. 

The pattern of genetic differentiation in Baetis 
mayflies has been studied within and among river 
systems (Monaghan et al. 2001, 2002, Hughes et 
al. 2003). Their studies were made in alpine envi- 
ronments that could be thought to give rise to geo- 
graphic isolation between drainages. In spite of 
this, the levels of genetic differentiation were in 
general low. This agrees with our results: allele 
freguencies were guite uniform within a species 
over geographic distances larger than the ones in 
the above studies. 

One may observe a horizontal pattern of allele 
frequency similarities among species at a single 
locality, as shown by the B. liebenauae sample 14 
and B. subalpinus sample 32 from site 11. We 
tried to get samples from the same sites in order to 
study sympatric populations and sampled the 
sites repeatedly. Nevertheless, B. macani and sp. 
nov. populations are not found to be sympatric. 
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4.2. What is Baetis macani Kimmins? 


Gill length is thought to be correlated with water 
flow in the environment, where the Baetis 
nymphs live. In stagnant water the size of gills is 
larger than in flowing water (Bauernfeind & 
Humpesch 2001). Baetis macani with broader 
gills and conspicuous tracheae lives in stagnant 
water (Savolainen & Saaristo 1981). Now it is ev- 
ident that this lentic form is a new species; here 
called sp. nov. The formal species description will 
be given in a separate publication. 

Our results show that north European B. 
macani (as recognized by Savolainen & Saaristo 
1981) with narrower gills consistently clusters 
separate from the forms with broader gills (Fig. 
2). The lotic form with narrower gills is identical 
with B. macani Kimmins since the nymph char- 
acters agree with the description of Macan 
(1957). We note that the nymph characters given 
by Lehmkuhl (1973) for the North American B. 
bundyae also agree with the description given by 
Macan (1957) for B. macani. Morihara & 
McCafferty (1979) pointed out that there are no 
consistent differences between B. macani and 
B. bundyae that would warrant a species status for 
B. bundyae. 


4.3. Conclusions 


Evidently more work is needed before we under- 
stand the taxonomy and ecology in this important 
group of mayflies. Phylogenetic methods tied to 
DNA analysis and F-statistics may solve much of 
the confusing overlap in morphological charac- 
ters and give us a good taxonomic understanding 
and good species designation that can be used in 
the ecological studies. Degerman er al. (1994) 
stressed that there certainly are difficulties in rec- 
ognizing the different forms when one is using 
Baetis in monitoring water quality. 
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